arXiv:1501.05033vl [math.NA] 21 Jan 2015 


NEW SIMULTANEOUS METHODS FOR FINDING ALL ZEROS OF A 

POLYNOMIAL 
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Abstract. The purpose of this paper is to present three new methods for finding all 
simple zeros of polynomials simultaneously. First, we give a new method for finding 
simultaneously all simple zeros of polynomials constructed by applying the Weierstrass 
method to the zero in the trapezoidal Newton’s method, and prove the convergence 
of the method. We also present two modified Newton’s methods combined with the 
derivative-free method, which are constructed by applying the derivative-free method 
to the zero in the trapezoidal Newton’s method and the midpoint Newton’s method, 
respectively. Finally, we give a numerical comparison between various simultaneous 
methods for finding zeros of a polynomial. 


I. Introduction 


With a typical iteration method such as Newton’s method, an initial approximation of a 
zero converges to a specific zero, but the Weierstrass method (or Durand-Kerner method) 
approximates all simple (real or complex) zeros of polynomial simultaneously (see |2, !]). 

Let P(z) = z n + ai£ n_1 + • • • + a n -\Z + a n be a polynomial of degree n having simple 
zeros with constants a±, a 2 ,..., a n . Let rr,..., r n be the distinct zeros of P(z) and let 
distinct complex numbers zi,...,z„ be their approximations. The Weierstrass method 
(Durand-Kerner method) is defined as 


„m +1 _ ro 


P(*r) 



(i) 
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for m > 0, and this method is one of the most frequently used iterative methods which 
give simultaneous computation of all zeros of P. If a function W t (z) is defined by 


Wi{z) 


P M 

n j¥i (z-zjY 


then W,j{z) has the same zeros as the polynomial P , and so the problem of finding the 
zeros of P reduces to that of zeros of the function W t (z). If we denote W t = W t (zi) for 
i — 1 , 2, ..., n in the case of z — © can be written as 


Bi = Zi — Wi , 


( 2 ) 


where z^ is a current approximation and z t is a new approximation to a zero of polyno¬ 
mial P(z). The method constructed by (J2J) is called the Weierstrass-like method (briefly, 
WLM). 

The aim of this paper is to present three new methods for finding all simple zeros 
of polynomials simultaneously. These new methods are based on the Frontini-Sormani’s 
midpoint Newton’s method ([7j) and the Weerakoon’s trapezoidal Newton’s method ([8]) 
which were modifications of the Newton’s method through iterative approximations. 

ft is well known that Newton’s method is defined by x* = x—jr^ with an approximation 
x and a new approximation x* of a zero, and is efficient to find a zero of an equation 
f(x) = 0 for a differentiable function / with proper conditions and a sufficiently close 
initial value (see m 

In [8j, Weerakoon proposed the trapezoidal Newton’s method defined by 

2 / 0 ) ^ 
X = X - - r---. Ml 

f'{x) + f\x*) 


He applied Newton’s method to the x* of the denominator. 

Along with ([2]), the midpoint Newton’s method that Frontini-Sormani proposed in [7] is 
constructed as 


x = x 


f(x) 


f (:X + \{x* - x)) 

They also applied Newton’s method to the x* of the denominator, and so set x* 


(4) 

f(P 
/'(*) ■ 
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Both the trapezoidal Newton’s method and the midpoint Newton’s method are of cubic 
order, while the original Newton’s method was of quadratic order. A variety of methods 
can be applied to the x* in addition to Newton’s method. Petkovic [5] derived the fol¬ 
lowing simultaneous method for finding all simple zeros of polynomials by applying the 
Weierstrass method to the x* in the midpoint Newton’s method: 


Zi = z t - 


P&) 


(5) 


P' {zi ~ \W t ) ’ 

which is called Newton- Weierstrss method (or NWM). Also, Petkovic (6j found the fol¬ 
lowing derivative-free method (or DFM) defined as 


i -p(z i -w i )/p(z,y (6) 

which has a similar form with the one above and this method is of cubic order. 

In this paper, we present three new methods for the simultaneous approximation of all 
simple zeros of polynomials by applying the Weierstrass-like method and the derivative- 
free method to x* in the trapezoidal Newton’s method and the midpoint Newton’s method. 

Throughout this paper, the convergence of zeros will be discussed and the order will 
be calculated for new constructed methods. We will use the notation a = Om(P) for two 
complex numbers a and b whose moduli are of the same order, that is, |a| = O (|6|). In 
addition, the error is defined as |e| = max {|e*|} with e i — z l — r t for i — 1 ,, n. 

2 — 1,...,72 

In all discussions, the order related to e,;, which is an error of the previously approx¬ 
imated zeros Zi, is presumed to be the same. After that, we will show that the order 
related to e*, which is an error of the approximated zeros concerning each method, is 
identical. For the same being, the order related to the already approximated zeros e \ is 
hypothesized to be identical as follows: 


e* = Om(g) for all i. 


In Section 2, we give a new method for hireling simultaneously all simple zeros of poly¬ 
nomials constructed by applying the Weierstrass method to the x* in the trapezoidal 
Newton’s method, and prove the convergence of the method. In Section 3, we present 
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two modified Newton’s methods combined with the derivative-free method. They are 
constructed by applying the derivative-free method to the x* in the trapezoidal Newton’s 
method and the midpoint Newton’s method, respectively. In Section 4, we give a numer¬ 
ical comparison between various simultaneous methods for finding zeros of a polynomial. 
Finally, we conclude that the convergence of all new constructed methods in this paper 
are similar or superior than other iterative methods of cubic order. 


II. Weierstrass-like Trapezoidal Newton’s method 


In this section, we construct a new method for finding simultaneously all simple zeros 
of polynomials of cubic order. By applying the Weierstrass method (J5]) to the x* in the 
trapezoidal Newton’s method (J3]) , we derive a new method constructed as follows: 


2 P( Zi ) 


(7) 


P'(zi) + P 1 {zi — Wj) ’ 


We call (J7J) the Weierstrass-like trapezoidal Newton’s method, and from this, simply, call 

it Method 1. 

The calculation and discussion of the order of Method 1 are similar to those of the 
Newton-Weierstrass method, which is an alteration of Petkovic’s midpoint Newton’s 
method (see 0)- From (0, we have the following theorem. 

Lemma 1. For a polynomial P(z), we have 



Proof. By the Taylor’s expansion around r*, we have that 




P" 
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From (ED, we obtain 

P(Zi) (P"(Zi)P( Zi 


P'(Zi 


2 P'{zi) 

P"( r i) 2 


2 0 «(e) + Om(s ) 


2V i+ 2P'(r i y i 


ef + Om(g 3 ] 


2/ l P "^C I I O ( C 3) X3 

1 P'( r <) ‘ 2 P'( n ) ‘ 1 } 


X 


P"{n) , P’"in) , P""(r*) 


p 


m + p%p + 2P^y e * + ° M(e } J ° M(e)+ ° M{e} 


— - (Om(c)) 2 (1 + Om{P)) 3 ^ P'(r-) ® M ( e ) ) 0«(e) + Om{& 3 ) — Chv/(e 3 ). 


□ 


From Lemma Q] we have the following theorem. 


Theorem 1. If the approximate zero ay grounded from Method 1 is close enough to ry 
and the order of e* is the same, then the order of ey is identical, and |ey| = Om{ |e| 3 ) 
formed. 


Proof. We easily see that the following equation is satisfied. 

n 

P( z i) = ]y[(^ - rf) = ( Zi - rf) - rf) = e* JJ(;y - rf) = 0 M {e). 

3 = 1 j¥=i I¥=i 


That is, 


w, = W,(z,) = O m (P(Zi)) = 0 M (e). 


( 9 ) 


If Q(z) = P(z) — n'=i(- — z j)j then Q(z) is a polynomial of order n — 1, and Q{zf) = 
P(zi) for all i. Therefore, Q(z) is the Lagrange interpolation of points Z\, z 2 , • • •, z n , and 
so we have 


j =i V kjtj j 

n 


z - z k 


X) h'illfe-^ IIP 


3 =1 





3 =1 
n 


Wi 


n 


Mi 


Zj - z k 




3=1 


k =1 
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Therefore, we obtain 


Wi 


Hz) = i + Etzt: n^-^) 


Z — Zj I 

3 =1 J / J=1 


fib 


w ’ n<* - z i) + 1 +1] ik* - z ^- 


j¥* 


... Z - Zj , . 

J / J=1 


From (gDD , it follows that 
P'(z) 


■d/* ( Z - Zj ) 2 


y _J_1 + 1 + Si* -(*-* ) E f w< 

m [jgz-ZjJ w, + (z-z i )(l + E^S- 


Substituting z = in (HID , we obtain 

P'(z t ) 


Wi 


P(zi) 


E l \ 1 s ^- Jr J+l Z .- Zj 




Zi - Zj 


Wi 


Therefore, we have 


W t = 


P(Zi) 1 + Zi - Zj _ P(Zi ) 

+ " ()) 


( 10 ) 


(ii) 


( 12 ) 


Now we will find the order of Method 1. If the Taylor’s expansion is applied to P'(zi — 
Wi), then we have 

2 P{zi) 

Zi Zi P'(zi) + P> (zi - Wi) 

_ 2 P( Z j) _ 

^ P'{ Zi ) + (P'(zi) - P"(zi)Wi + \P'"{zi)W? + • • •) 

= ,_ EM _ 

* P'(zi) - 1 P"{zi)W z + \P'"{zi)Wf + • • • ■ 
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By OH]) , (TT2j) and Lemma [TJ we have that 

P&) 


Zi = Zi 


P'{zi) ( X -S^Wi + OMie 2 )) 


= Zi - 


= Zi 


1 + 

1 + 


P"(Zi 


2P f {z., 

P"(Zi)P{Zi) 

2P'( Zi ) 2 


= Zi — 


= Zi 


P(z) 

P'&) 

P(Zj) 

P'{Zi) 

P{Zi) P"(Zi)P(Zi ) 2 P( Zi ) (P"(Zi)P(Zi) 


Wi + 0 M (e 2 ), 

(1 + Ojw(e)) + Om{& 2 , 


P'(zi ) 2 P'{ Zi f P '{ Zi ) V 2 P'(ziY 


-Om(g) + Om(c 2 ) 


p(zi) p"^^ 


P'{zi) 2 P'{ziY 

The Chebyshev’s method is defined by 

x = x — ( 1 + 


Om(s 3 )- 


f'(x)f(x)\ f(x) 


2 f'(x) 2 J fix)’ 

and of cubic order (see [7] Section 5.2]). According to the Chebyshev’s method, we see 


that 


a - -S4 - ALAM! - r < = 0„(e 3 ). 


p\z,) 2 p\z,y 

Therefore, the order of e z is calculated as follows: 


eh = Zi - n = 


*• P'(z t ) 2 P'(z,y P + °M(e) 


— Om (l e | 3 ) 


□ 


III. Modified Newton’s methods combined with Derivative-free method 


In this section, we present two modified Newton’s methods combined with the derivative- 
free method flSD for finding all simple zeros of a polynomials simultaneously. The one is 
a form that the derivative-free method is applied to the x* in the trapezoidal Newton’s 
method (J2J) as follows: 


Zi = Zi 


2 P(zi) 


P'(Zi) + P' Izi- — 


Wi 


-P{zi-Wi)/P{zi) 


(13) 


which is called the Derivative-free trapezoidal Newton’s method, or simply, Method 2. 
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From (TT3]h we have the following theorem. 

Theorem 2. If the approximate zero x t grounded from Method 2 is close enough to r l 
and the order of e* is the same, then the order of e^ is identical, and |e*| = Om( |e| 3 ) is 
formed. 


Proof. Since Petkovic’s derivative-free method (J6]) is of cubic order 

W, 3 . 

n = 0 M (e ), 


(14) 


1 - P(zi - Wi)/P(zi) 

(see ®)- Using (dU) and the Taylor’s expansion, £$ is calculated as follows. (In this case, 
r< _ 1 

“ j\ ' P'(n) >■ 

2 P{zi) 


Zr = Zi 


= Zi 


= Zi 


P'{Zi ) + P' (zi — — 


Wi 


P(zi-Wi)/P(zi) 


2 P(ri + e-i) 


= Zi 


P'(ri + ef) + P' (n + 0 M (e 3 )) 

2 (P(n) + P\n) ei + ±P"(r, : )e 2 + \P"\ri)e 3 + • • •) 

(Pfn) + P"{n)ei + \P"\rf)e\ + ••■) + (P'(u) + P"(r*)0 M (e 3 ) + ■ • ■) 
2P'(ri) (e* + C* 2 e 2 + C^ef + Om(c 4 )) 


2 P'{rf) (l + C 2 e* + |C 3 e 2 + Om(g 3 )) 

— Zi — e* (l + C' 2 e,: + C^e 2 + Om(c 3 )) ^1 — C’ 2 ej + -C^e 2 + ^M(e 3 )j 

= ^ - e* + Om(c 3 ). 


Therefore, the order of e, : is calculated as follows. 


= - u| = |z* ~ei + 0 M (e 3 ) - Vi | = 0 M (|e| 3 ) . 


□ 


Now we apply the derivative-free method to the x* in the midpoint Newton’s method 
d3J) and construct the iteration as follows: 


% = Zi - 


P(Zi) 


P' Zi-h 


Wi 


2 1 -P(zi-Wi)/P{zi) 


(15) 
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which is called the Derivative-free midpoint Newton’s method. From this, we call it 
Method 3 simply. From (fT5]h we have the following theorem. 

Theorem 3. If the approximate zero Xi grounded from Method 3 is close enough to ri 
and the order of ej is the same, then the order of e* is identical, and |e*| = O m (| e| 3 ) is 
formed. 


Proof. By using and Taylor’s expansion, i* is calculated as follows: (In this case, 
r> — 1 P u Hn) \ 

W “ j\ ' P'( n ) >■ 

* _ ._ pm _ 

P' ( z . — I . _ Wi _ ) 

2 1 -P(zi-Wi)/P(Zi)) 

= z _ n*) 

P' (■\zi + \p + 0 M (e 3 )) 

_ P{n + ej ) 

P' (^ + \ei + 0 M (e 3 )) 

P(ri) + P'{ri)ei + ^"{r^ej + \P"'{ri)e f H- 

P'(rf) + P"{rf) Qe* + 0 M (e 3 )) + Qe* + 0 M (e 3 )) 2 H- 

_ PfPi) (Cj + + C3 e i + Om(c 4 )) 

“ * ~ P'fc) (1 + + fc 3 e? + 0 M (e 3 )) 

= Zi — ei (l + C 2 ei + C^ef + Om(c 3 )) ~~ C 2 e, h + — + Om(c 3 )^ 

= 2* - e* + 0 M (e 3 ). 

Therefore, the order of e* can be calculated as follows. 


ei | = |£j - ri| = \zi -ei + 0 M (e 3 ) - rf\ =O m (|e| 3 ) . 


□ 


IV. Numerical comparison 

In this section, we give numerical experiments and comparisons between various simul¬ 
taneous methods for finding zeros of a polynomial. These methods are all of cubic order. 
They include Method 1, Method 2, Method 3, the derivative-free method (DFM), the 
Petkovic’s Newton-Weierstrass method (NWM), and the Weierstrass-likc method (WLM). 
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For a polynomial P(z) = z n + aiz n * + • • - + a n -iz+a n , we choose initial approximations 
as Aberth’s approach (see my 

(I (“-!))■ 


In this case, R is a radius of a circle, where the initial zeros by Aberth’s approach are 
located in complex number plane. We use the following Henrici’s fomula to select R (see 

w- 


R = 2 max |aJ 1/fc . 

l<fc<n 

According to Henrici’s formula, a disk {z : \z\ < R} centered at the origin contains all 
zeros of polynomial P(z). 

The polynomials that we used on numerical comparison are as follows: 


Pi — {x — l)(x — 2)(x — 3)(x — 4) 

P 2 = (x — l)(x — 2)(x — 3)(x — 4)(x — 5) (16) 

P 3 = (x — l)(x — 2) (x — 3)(x — 4) (x — 5)(x — 6) 

P 4 — x 8 + 5x 7 + 3x 6 + 7x 5 + 6 x 4 + 8 x 3 + x 2 + 3x + 7 


Here Pi(x), P 2 (x) and P 3 (x) are Wilkinson’s polynomials when n — 4, 5, 6, respectively. 
We approximated the zeros until it satisfy the following condition: 


max 

l<k<n 


P(z, 


W' 

k / 


< 10 


-10 


(17) 


In Table I, we give a numerical comparison between several methods to find all zeros of 
those polynomials (116]) . It contains the iteration number m and the value max P(z^ n> ) 
of iterative methods, after we approximated (TTT)) to a satisfying label. The smaller the 
m, the faster approximated on the zeros, When m is the same, it can be interpreted that 


a smaller max 

l<k<n 


p\~, 


(™)' 


leads to a higher accuracy of approximation. All computations 


have been done using MATLAB. 
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Table I. The number of iterations (the error) of iterative methods 


Poly. 

(USD 

Method 1 

(ED 

Method 2 

03 

Method 3 

(ESD 

DFM 

m 

NWM 

© 

WLM 

© 

Pi 

9(8e-14) 

8(9e-14) 

7(le-14) 

9(9e-14) 

8(le-10) 

13(3e-12) 

P 2 

12(7e-13) 

ll(4e-13) 

9(4e-13) 

ll(2e-13) 

ll(le-12) 

17(2e-12) 

P 3 

14(2e-ll) 

13(5e-12) 

ll(8e-12) 

13(8e-ll) 

13(7e-12) 

21(2e-ll) 

Pa 

14(2e-ll) 

13(2e-ll) 

10(2e-ll) 

14(2e-ll) 

13(2e-ll) 

21(2e-ll) 


V. Conclusion 


In this paper, three new methods for the simultaneous approximation of all simple 
zeros of polynomials by utilizing the trapezoidal Newton’s method and the midpoint 
Newton’s method were proposed. It was proven that each method was of third order. By 
simultaneously approximating all simple zeros of polynomials and by comparing numerical 
experiments with various methods that are of third order, we obtained that the results 
of Method 1 and Method 2 are similar with that of previous methods. But we found out 
that the result of Method 3 are superior than that of any other methods. All methods 
we constructed in this paper are new and creative. It seems that these methods can be 
applied to various fields, and the study on the applications of Method 3 is now in progress. 
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